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Physics of Systems Containing Neutron Stars 


NASA Grant NAGW-1618 


Jacob Shaham, the PI on this grant, died tragically at the end of April 1995. He was 
also the sole investigator and his achievements on it are not easily summarized by anyone 
not involved in his work and familiar with his final thoughts about it. The sudden onset 
of his incapacitating illness began shortly after his submission of his last (January 1995) 
progress report. This final report is based mainly upon a reading of it and access to his 
3-year renewal proposal. As much as possible it incorporates his own words. 

* * * 

This grant dealt with several topics related to the dynamics of systems containing a 
compact object. Most the research dealt with systems containing Neutron Stars (NSs), but 
a Black Hole (BH) or a White Dwarf (WD) in situations relevant to NS systems were also 
addressed. Among the systems were isolated regular pulsars, Millisecond Pulsars (MSPs) 
that axe either Single (SMPs) or in a binary (BMPs), Low Mass X-Ray Binaries (LMXBs) 
and Cataclysmic Variables (CVs). Also dealt with was one aspect of NS structure, namely 
NS superfluidity. 

A large fraction of the research dealt with irradiation-driven winds from companions 
which turned out to be of importance in the evolution of LMXBs and MSPs, be they SMPs 
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or BMPs. While their role during LMXB evolution (i.e. during the accretion phase) is 
not yet clear, they may play an important role in turning BMPs into SMPs and also in 
bringing about the formation of planets around MSPs. 

Work was concentrated on the following four problems: 

The Windy Pulsar B 1957+20 and its Evolution 

In PSR B 1957+20, the observed orbital period V seems to change over timescales T 
of order 10 7 — 10 8 yrs (Arzoumanian et al., 1994). If £ is the angular momentum per unit 
mass for the binary and if the total gravitational potential is <£>, for a small companion 
mass and for a small deviation from ordinary point-mass gravity the change in total energy 
8E can be written as 

6E = e 0 Sm-F^-+H^, ( 1 ) 

where H is a constant, e Q measures the excess specific energy in the wind and where, if 
the part of $ that represents deviation from point masses is roughly a power law in a, 
F ~ 0(m<I>) = 0(mv or b 2 )- It follows from Eq. (1) that the corresponding change in 
orbital velocity of the binary over the above timescale involves a rate of energy change 
(lost or gained) of ~ mwort+orb ~ Grn J£vs ~ io 31 — 10 32 erg/sec. If a wind is the broker 
for these energy exchanges (note that the Roche lobe [RL] surface in B1957+20 intercepts 
6 x 10 32 erg/sec from the pulsar), then m ~ 10 16 — 10 17 g/sec (since Uorb 2 ~ 10 15 cm 2 /sec 2 ). 

By combining self-consistently the idea of magnetic activity on the companion of 
B 1957+20 (Applegate 1992) with the idea of the intense irradiation- activated mass loss, 
Shaham and Applegate put together a scenario for the system that seems to be able to 
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account for many features of B 1957+20 (Applegate and Shaham 1994). In this scenario, 
magnetic activity directs the wind to leave the surface on field lines essentially out to 
the eclipse shock front. As a result, the companion spin is being magnetically breaked 
continuously. Tidal torques from the NS act to restore corotation and deposit entropy in 
the companion via tidal friction. This heat is sufficient to support the companion at its 
bloated configuration as well as to power its “dark”, side luminosity in steady state, thus 
making the companion of B1957+20 the first discovered tidally-powered star. 

From the above it can be determined when the evaporation in B1957+20 will be 
expected to really take off to very high mass-flow rates. If the companion continues to 
fill the same fraction of its RL, then further evolution will cause, for some companion 
mass and at the present pulsar luminosity, the surface temperature to exceed the escape 
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temperature. That that mass is: ~ 40 M®a p 4 , with a p the orbital separation in units of 
the present one. It is easy to speculate, by analogy with Jupiter, that magnetic activity 
will still go on at that stage so that the companion will clearly maintain a fixed ratio with 
its RL. If it does not, we would simply have to wait for a lower mass to have the onset of 
intense evaporation. Once it has begun, however, the companion will evaporate completely 
because the illumination will provide super-escape thermal velocities throughout the whole 
process. 


Wind “Echoes” in Tight Binaries 

Several of the Soft X-Ray Transients (SXTs) from BH candidates were seen to be 
followed by secondary outbursts. A good example is Nova Persei, GRO J0422+32, that 
erupted in August of 1992 and by the end of 1993 had four secondary outbursts; the earlier 
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ones were seen by BATSE and ASCA, some of these and the later ones were looked for 
and seen in the optical. The overall decay of the x-ray luminosity was typical: an almost 
perfect exponential decay with an e-folding time of order of a month, interrupted by the 
secondary bursts at roughly 4 month intervals. The optical bursts were less clear in their 
luminosity trend even though they did keep the inter-burst intervals and roughly coincided 
with the x-ray maxima when x-ray data was taken; however, in all bursts except, possibly, 
for the last one (of Dec ’93), the optical luminosity may have been but a small fraction of 
the bolometric luminosity. 

In modeling the SXT phenomenon one may want to separate the main burst from the 
secondary ones. Chen et al. (1993) pointed out that the latter might reflect some response 
of the system to the main burst. Following the first two secondary bursts, Shaham and 
collaborators (Augusteijn et al. 1993) suggested a simple model for the whole pattern of 
secondary bursts in terms of successive wind “echoes” of the main burst, which was able to 
predict correctly the times of the last two secondary bursts and the bolometric luminosity 
of the first of those (it was not possible to look at Nova Persei with a sensitive-enough 
x-ray satellite during the last secondary burst). It is in the framework of this model that 
they continued theoretical research on wind “echoes”. 

In the model, they assumed that the centrally emitted x-rays are able to illuminate the 
companion through scatterings in a large corona, without having to worry about obscura- 
tions in the disk itself or in an inner corona. As the main burst occurs, it induces a burst of 
extra mass flow from the companion. This extra mass goes first into an orbit which is not 
too different from that of the companion, because the ratio of companion-to-BH masses is 
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quite low. It takes several months for the material to drift into the inner disk, where x-ray 
emission occurs, and several days to drift through that region into the BH. The new burst 
of x-ray emission has less total energy than the main burst, but it nevertheless lifts a new 
burst of mass flow from the companion, that begins its way again into the BH to produce 
the next x-ray burst. All of the above assumptions guarantee that the secondary burst 
pattern comes about by simple linear response in the system with constant parameters. 

As they predicted, J0422+32 returned back to its quiescent state. The optical data 
collected during the “echoes” suggest that the disk is getting cooler as the overall intensity 
(hence rh?) drops, since the optical luminosity does not seem to drop from one “echo” 
to the other as fast as the X-rays did during the first “echoes” (when X-Rays were still 
visible). A crucial part of the understanding of the system (and the several others that 
seem to have been detected before and since Nova Persei) will be in modeling the “echo” 
to “echo” variations in luminosities of the various energy bands, while separating contri- 
butions from reprocessing of radiation in higher energy bands from direct disk luminosity, 
and in modelling the corona. Recent calculations (Murray et al., 1994) indeed show that 
multiply-scattered radiation does act to partially offset the loss of direct radiation from 
the central source, predicted by earlier models to occur because of the inner corona; these 
multiple scatterings could well bring the radiation to the companion regardless of any 
close-to-the-plane shadowing effects. 

In their phenomenological model (Augusteijn et al., 1993), the important physics was 
contained in the interburst time T and in the disk spread function: they assumed that 
after spending the time T at a distance of order of the companion distance matter goes 


5 



through the inner disk on a timescale of order /3 — 1 , while the flow stretches over a similar 
time scale (see, e.g., Bath et al., 1974). As they found from their fits to several bursts 
(notably GS200+25, 0620-003 and Nova Persei), /3 -1 <C 7 -1 , where -y~ 1 is the main burst 
decay time; hence the detailed shape of the spread function is, actually, not very critical 
for this picture. 

In the original suggestion Shaham et al. we did not comment on why the decay of 
the main burst was so clearly exponential, as that did not have anything to do with the 
echoes. They have begun to look into this question. 

Over the years, there have been two general frameworks in which the main burst was 
modeled: either as an inner disk instability (Tuchman et al., 1990) or as a companion 
mass transfer instability due to accumulation of excess entropy in the companion during 
the quiescent state (Hameury et al., 1986). The good fits that they obtained for SXTs 
with our simple linear model described above are actually suggestive that the second of 
these options is the right one. One should understand how the exponential decay of the 
main burst comes about, assuming the latter is, indeed, due to a mass flow instability. 

It is interesting that the e- folding time for the exponential decay of the main burst is 
of the same order (~ 1 month) as the interburst interval (~ 4 months). If the interburst 
interval represents, indeed, the time it takes matter to get from the LI point down to the 
BH, it seems natural to interpret the e-folding time as due to matter motions inside the 
RL of the BH as well. Shaham et al. argued, quite generally, that in a stably accreting 
system, matter torques (coming from the mass flow stream and the disk, say) control the 
stability of the steady-state value of m. Thus, if m drops suddenly, the change in tidal 
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torques due to the new flow parameters will be such as to torque the binary into getting 
closer, so as to make the companion fill its RL a little more. Inversely, when a sudden rise 
in ra occurs, the extra flow will change the tidal torques so as to make the binary expand 
slightly and make the companion move further from RL contact. During a transient event, 
the excess mass flow rate, m e , will therefore depend on the excess torque, hence (while 
excess matter still accumulates in the BH RL) on the total accumulation of the excess 
mass flow, m e dt (we assume that the event begins at t — 0). This reasoning leads to an 
equation of the type 


rh e (t ) oc (-) 



(2) 


hence the exponential behaviour. 

Shaham had begun work to find Eq. (2) in simulations of the flow during non-steady- 
state situations. 


Post Nova X-ray Emission in CVs 

ROSAT observations of Nova Mus ’83 (GQ Mus) showed some ten years of strong 
soft x-ray emission (Ogelman et al., 1993; x-rays seem to have turned off only recently). 
The spectral characteristics, as well as the Eddington-magnitude luminosity, were quite 
reminiscent of the LMC sources CAL 93, CAL 87 and RXJ 0527.8-6954, which are thought 
(van den Heuvel et al., 1992) to contain WDs on which nuclear burning of accreted matter 
occurs. Stable nuclear burning of this type was found to occur for mass accretion rates 
between 10 7 and 4 x 10 7 MQyr 1 , which can be driven by either GR or Magnetic 
Breaking (MB) off 1.4 — 2.2 M© companions. The orbital periods in question here being 
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1.04 and 0.44 days respectively for CAL 83 and CAL 87, this constitutes a consistent 
picture for these sources. 

In spite of the similarities, GQ Mus could not possibly fit into the above picture, 
because its 85.5 min orbital period implies a ~ .1 Mq companion (if the companion is to 
be close to filling its RL). Shaham suggested that a decade-long episode of “self-excited” 
mass accretion could have taken place here. While the feedback factor here is, apriori, 
around 10 3 less than it would be were the compact object an accreting NS and were the 
x-ray luminosity due to free fall energy, it is only a factor of ~ 30 smaller if the x-rays are 
produced by nuclear burning of accreted material. 

There is some evidence that accretion was present in the system while the x-rays 
were on. Optical observations reveal (Diaz and Steiner 1990) emission-line velocities and 
a light curve that is consistent with an AM type binary, whose companion is illuminated 
by the WD. For a companion radius R and surface temperature T, R 2 T is constrained 
by these observations to be around 6 x 10 24 cm 2 AT, yielding a surface illumination of 

4.5 x 10 35 jRio ~ 6 erg /sec, with i?i 0 = f o ™ cm (for a main sequence companion of .16Mq, 
Rio 6 ~ 2). From the limit on the solid angle by which the WD sees the companion, the WD 
luminosity can be placed at > 3 x 10 37 erg/ sec, thus confirming the x-ray estimates (which 
are always somewhat uncertain because of the high extinction factors at these wavelengths). 
Shaham and collaborators modeled the Diaz and Steiner light curves in more detail and 
found them consistent with a RL-filling companion that is strongly illuminated by the WD 
but shadowed by an AM Her type stream of mass flow from the LI point. 

If bootstrapping really happens in the post Nova GQ MUs, the efficiency of illumination- 
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to-mass-flow conversion would be similar to that needed to explain BMP B1957+20, an- 
other source with very high such efficiency. Shaham had begun exploring the possibilities 
for this particular system and for post-Nova CVs in general. 


Dynamics of Pinned Superfluids in Neutron Stars 


Neutral superfluids can rotate if they contain quantized vortices of density 


2 n 


n ' V = T~7'10 6 Vl 0 3 sec _ 1 


ft 


cm 


-2 


( 3 ) 


where k is the quantum of vorticity, k ~ 3 • 10~ 4 cm 2 sec _1 . Charged superfluids (i.e. 
superconductors) can carry a magnetic field if they contain quantized magnetic vortices of 
density 




( 4 ) 


where <f* 0 ~ 2 - 10~ 7 Gem 2 is the flux quantum. 

Quantum fluids form in NS interiors by Cooper pairing of “dressed” crust or core 
neutrons or “dressed” core protons (Sauls 1989). Cooper-pair binding energies are a few 
hundred keV and depend on the local density and composition of matter. Cooper pairing 
is broken in the vortex cores, which have to contain normal fluid; therefore vortices will 
form locally where binding is weakest to begin with. This will create an affinity (pinning) 
of vortices to certain microscopic locales in the NS such as crustal nuclei, spaces between 
nuclei, or cores of vortices of the other kind (in the NS core). Pinning forces are of order 
of the binding energy gradients, ~ 10 17 — 10 dyne /cm /vortex. 

The global rotation of a superfluid is determined by the locations of its vortex lines. 
These, in turn, move downstream unless a force acts on their normal cores; in this case they 
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acquire an extra velocity component perpendicular to the force, of a magnitude proportional 
to it. When embedded in a normal fluid that rotates slower than the vortex density implies 
[by Eq. (3)], tangential friction between the normal fluid and the vortex cores sends the 
vortices outwards and the vortex lattice dilutes itself for the lower rotation rate. However, 
when pinned to centers corotating with the normal fluid, the vortices will not move unless 
the Magnus force overcomes the pinning. In general, friction will produce a microscopic 
displacement of the vortex radially outwards, into the pinning energy gradient region. The 
resulting radially inwards pinning force will introduce a backward tangential velocity for 
the vortex , thus forcing the vortex to move with the (lesser) speed of the pinning centers, 
in spite of the higher vortex density. The superfluid will only be able to slow down if 
the Magnus force (proprtional to the velocity difference between the vortex core and the 
superfluid) causes unpinning or if it can break the lattice and carry the pinning sites with 
it. 


Sudden unpinning events have long been considered good candidates for causing spin 
up “glitches” in pulsars (see, e.g., Alpar and Pines 1989): the sudden slowdown of the 
superfluid is countered by the observed sudden spinup of the crust. One of the mysteries 
of this scenario is how a collective unpinning event comes about. One starts with a single 
vortex, unpinning under the appropriate maximal Magnus force. As it is spiralling out, a 
superposition of the motion downstream and the radial component introduced by friction, 
it collides with other pinned vortices which are also on the verge of being unpinned. The 
extra (small) velocity that it induces on them during this close encounter can bring about 
their final unpinning and so, in an avalanche, a bunch of unpinned vortices is formed. 
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Another mystery has to do with the distance which a bunch travels until repinned, if 
repinned. To repin, a vortex must lose its energy to something. Back-of-the-envelope 
calculations show that the energy of an unpinned vortex under typical conditions fits right 
into the energy gap between acoustical and optical lattice phonons and it would not be able 
to lose its energy to the lattice. Shaham considered energy loss of the vortex to friction as 
the main repinning mechanism. Both the travel distance and the initial bunching determine 
the observed rise time and magnitude of a “glitch”. 

The basic equation governing the motion of a single vortex line, in the rigid vortex 
approximation, is 

p s [(v s - v L ) x k] - tj(v l - v n ) + VfK ~ v n ) x k] + f P = 0 , (5) 

where p s is the superfluid density, rj and rf the longitudinal and transversal viscosity 
coefficients, k the vortex vorticity (of magnitude k), v s the local superfluid velocity, 
the velocity of the vortex line, v n the velocity of the normal component and fp the pinning 
force. All bold-faced quantities are 2-D vectors in the plane perpendicular to the rotation 
axis except for k, which is parallel to that axis. In the above equation, the first term is the 
Magnus force and the next two terms represent the longitudinal and transverse frictional 
forces due to the interaction with the electrons and the crustal phonons in the crust. The 
pinning force fp is a function of the relative vector between the vortex line and the pinning 
site: if r is the location of the vortex and R z the location of the i — th pinning site (assumed 
here to really be a pinning line), then fp = fp(r — R t ). 

A second set of equations describes the motion of the pinning sites. These are acted 
upon by the ordinary forces in the NS crust as well as by the pinning force due to the 
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vortex and can be simply represented as a coupled set of harmonic oscillator equations or, 
for more simplicity, by a single harmonic oscillator equation with some angular frequency 

CJq . 

At the time of his death Shaham was working on the assumption that v n is given and 
that, as a first approximation, so is v^. At a later stage he intended to put in by hand some 
“glitch” function to mimic self consistently the superfluid rotation based on the picture 
based upon the behaviour of a single vortex. The binding energy would be first modeled 
as some Fermi function which vanishes at the center of the pinning site and rises to the 
value of ec,max ~ £c,min, the maximum variance of Cooper pair binding energy per particle 
values ec, over a distance of order 10/m (the force is the gradient of this function). It was 
then planned to calculate it more accurately from the most recent local equations-of-state 
around crustal nuclei. 

Unpinning occurs when the energy density of relative motion grows to a sufficiently 
high value, i.e. to roughly €c,max ~ £c,min in the region. Lattice breaking occurs under 
similar conditions involving the lattice. This limiting energy density may well exceed in 
some regions the value of ec.max, the binding energy per particle at the pinning site. In 
this case, it may be energetically favorable to form a new vortex of the opposite (negative) 
vorticity on a neighboring pinning site and decrease the positive vorticity in this way 
rather than move (the positive) vortices out. Vortex formation is a local process. Since a 
mismatch exists between the motions of the superfluid and the (slower) normal component, 
the latter is seen locally to rotate in the reverse direction; so, with the availability of energy, 
a negative vortex is sure to form if vortices can at all form in the region. 
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If negative vortices do form, many ideas regarding NS crustal superfluidity will have to 
be reexamined. These include unpinning events as sources of “glitches” , and magnetic field 
decay during crustal “continental drift” induced by pinned vortices breaking the lattice 
and moving outward carrying the lattice with them (Ruderman 1991). Shaham began to 
assess the possibility of forming negative vorticity in NS crusts in detail. To determine the 
possibility for negative vortex formation he needed good estimates for Cooper pair binding 
energies in the local environment of a crustal nucleus. While these do exist (and suggest 
favorable conditions at least in some parts of the crust), they may not be reliable enough 
to detemine whether the energetics favors negative vortex formation. It is hard to see how 
to do better calculations on this question, because the main uncertainty are the many- 
body calculations, which were not tested in the laboratory for these systems. Secondly, he 
needed to determine how these vortices form even if the energetics is favorable, because 
vortices can only form in regions bordering normal neutron fluid. 

One way to form them is this: an existing vortex may stretch and elongate in a non- 
pinned segment, first up the axis to be equivalent to a negative vortex and then outwards 
and down the axis again. This will effectively create a negative vortex next to the pinned 
existing one, as well as another positive vortex further away. One needs to create two 
normal core lines here, not one (of course they should end up pinned), so the energetics 
will be different, but if repinning is prevented until the full new lines are in place, it will 
be a viable process. Alternatively, repinning may take place continuously and each gener- 
ation would put out additional loops; this would take more than two new normal cores per 
existing vortex but it is only the instantenous energetics that is of concern: there is enough 
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energy in the rotational energy reservoir to form many vortex cores per each existing one. 

An MSP that was spun up from a low spin while its interior was already at tempera- 
tures much below the superfluid transition is going to need to fill itself up with a number 
of neutron vortices that far exceeds the number it started with, to accommodate to the 
new spin rate. Otherwise its superfluid neutrons are going to continue to rotate at their 
old rotation rate. Any way which the superfluid can find to form these vortices will also 
be a way for it to form negative vortices once slowdown begins (provided the enrgetics 
is right). Thus, either MSPs do not have enough neutron vortices to have any influence 
on the magnetic field or they do, in which case they may again not have any effect on B 
during spin-down because negative vorticity will form. 
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